Bosonic model with Z3 fractionalization 
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Bosonic model with unfrustrated hopping and short-range repulsive interaction is constructed 
that realizes fractionalized insulator phase in two dimensions and in zero magnetic field. Such 
phase is characterized as having gapped charged excitations that carry fractional electrical charge 
1/3 and also gapped Zi vortices above the topologically ordered ground state. 



I. INTRODUCTION 

A flurry of recent theoretical activity has produced spe- 
cific model system pftalizations of fractionalized phases in 
two dimensions flla'afflu Essentially all of tlie fractional- 
ized states constructed so far are states.Q On a formal 
level, these realizations employ the following route to Z2 
fractionalization: Strong local correlations lead to a C/ (1) 
gauge theory as a low-energy description; this gauge the- 
ory is then driven into a deconfined state by a condensa- 
tion of objects carrying gauge charge 2 that also appear 
in the low-energy description. This formal structure has 
been brought out very directly in Refs. |5[^,§. On a more 
physical level, the fractionalized insulator is produced de- 
parting from a supieEGpnducting state by a condensation 
of double vortices.Bliy The main body of work concen- 
trated on the Z2 states since these are expected to be 
the simplest to realize. However, it is clear that more 
complicated fractionalized states are also possible. For 
example, it is conceivable that in some system the super- 
conducting state is quantum-disordered by a condensa- 
tion of triple vortices; the resulting insulator is then a Z^ 
fractionalized state. 

In this paper, we indicate how a Z3 fractionalized state 
can be engineered in a relatively simple bosonic model 
with unfrustrated nearest-neighbor hopping and short- 
range two-body repulsive interaction. Much of the con- 
struction parallels closely the Z2 examples of Refs. ||,|[^: 
The low-energy Hilbert space is selected — by stipulating 
particular charge interactions — in a manner that natu- 
rally admits splitting boson charge into three pieces; this 
Hilbert space is protected by a large charge gap. The ef- 
fective description of the fractionalized state has gapped 
chargons carrying electrical charge 1/3 and coupled to 
some special Z^ gauge theory which we analyze in detail. 
Our main message here is that one does not need very 
contrived systems to obtain more complicated fraction- 
alization patterns. 



II. Z3 VIA CHARGE FRUSTRATION 



FIG. 1: Josephson junction array formed by penetrating 
hexagonal (r) and triangular (i?) lattices, modeled by the 
Hamiltonian Eq. (^). The shaded area indicates schemati- 
cally hexagon charging energy UN^. 



bosons residing at the hexagon centers (upper-case R). 
Bosons can hop between the neighboring sites as indi- 
cated by the links on the figure; the hopping amplitudes 
are wi for (r, R) links and W2 for (r, r') links. We also 
stipulate strong repulsive interactions that favor charge 
neutrality of the hexagons, in addition to the on-site re- 
pulsion that favors charge neutrality of the individual 
sites. The complete quantum rotor Hamiltonian is 



R,reR 



{rr') 



Here, 



The model is defined on the lattice shown in Fig. 
which we can think of as a hexagonal lattice with ad- 
ditional sites placed at the hexagon centers. We have 
?/)t = e.^i'r bosons residing on the hexagonal lattice (these 
sites are always labelled lower-case r), and h^p, = e*^^ 



^r} Eire conjugate number-phase variables 
[e.g., in the phase representation n\ = —id/d9ii, bl^tpr + 
h.c. = 2cos(6'fl — (pr)], and similarly for {nf,(j)r}. The 
number-phase variables are particularly appropriate if we 
think of the model as describing an array of Josephson- 
coupled superconducting islands. Both bn and ipr bosons 
carry electrical charge qb- 

In the above Hamiltonian, Nn is the boson number 
associated with each hexagon: 



R 



re-R 



,4- 



(2) 
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Thus, the total boson number in the systenJli is 

General analysis of the possible phases proceeds as 
m Ref. ||. Here, we focus on the insulating states 
that are obtained for large U 3> wi,W2,Ub,u,p. If the 
'Wi,W2,Ub,U',l, terms are all zero, there is a degenerate 
manifold of ground |States specified by the requirement 
Nr = for each The ground state sector is sepa- 
rated by a large charge gap U from the nearest sectors. 
Including the wi,W2,Uh,u^ terms lifts this degeneracy 
in each sector, which is best described by deriving the 
corresponding effective Hamiltonian in each sector. 

The effective Hamiltonian in the ground state sector 
{Nn — 0) is, to third order in the perturbing terms, 

r 

-^ring 51(^1 V'2V'|V'4V'JV'6+ll.C.) . (4) 

o 

Here, stands for the on-site repulsion terms as 

in Eq. ([I]); i?i,i?2,i?3 label the three hexagon centers 
adjacent to r; = wJ/(6C/^) and i^Hng = 3w|/[/^. 

This is our main step in obtaining fractionalization. 
The above Hamiltonian looks similar to a compact U{1) 
gauge theory coupled to a charge 3 scalar field. Thus, if 
we think of the ipr as some gauge fields, then it is very 
suggestive to think of the as carrying gauge chargej-S. 
[see also Eq. (^]. Standard Fradkin-Shenker analysisE^ 
then suggests that by condensing the bn field, which can 
be arranged by making Jc, i^ring large, we can deconfine 
objects carrying gauge charge 1. By virtue of Eq. (||), 
such objects carry fractional electrical charge (7b/3, and 
we obtain a fractionalized insulator. 

The reader who finds the above statements believable 
may now declare victory in achieving fractionaliza- 
tion. However, if we want to describe the deconfined 
phase(s) in any detail, we need to study the above Hamil- 
tonian directly since it is not related in any simple man- 
ner to the conventional gauge theory. This is our focus 
in the remainder of the paper. Of course, we will confirm 
the deconfinement, but we will also find that the decon- 
fined state that obtains for large if ring on all hexagons is 
in fact & Z^x Z3 state (see below for details). 

Proceeding with this analysis, consider the regime of 
large Jc ^ 00 and small u;, — > 0. It is convenient to 
perform the following change of variables. Define the 



operators h\j^ = e" 



and V^J 



'cR2 0cR3 



(5) 



Here, sr^I, e'^''/^, or e"^''/^, so that OcR e [0, 2tt). One 
can readily verify that Nr is conjugate to 9cr, while nf 
is conjugate to cpr- carries electrical charge q^/i and 
can be thought of as a chargon field, while V'J is charge 



neutral. These new variables are indeed natural in the 
description of the deconfined phase, but to recover the 
physical Hilbert space, we need to impose the constraint 



exp 



27r 



r£R 



(6) 



In the new variables, the Hamiltonian becomes 

r r 
-■faring E ('^lV'2V'3V'4'05V'6 + h- 



O 



(7) 



Note that we are left with gauge fields ipr only, since we 
are working in the "uncharged" ground state sector 

The term Jc acts as a Z3 anisotropy on the c/fr field. In 
the limit Jc 00, e"^'' becomes a Z^ field: (pr — 0, 27r/3, 
or 47r/3. The operator P+ = ezli^-^/^X shifts the states 
of the quantum Z3 clock by -t-l,L3 whereas the constraints 
specifying the uncharged sector Nr = become 



= 1 



(8) 



r£R. 



The u^{nf)'^ term causes tunneling between the different 
states of the quantum clock, and this can be described 
by an effective "transverse field" —h{P^ + Pr) acting on 
the clock. 

The effective Hamiltonian now becomes Z3 "ring- 
exchange" Hamiltonian on the hexagonal lattice of r sites 



ifri 



O 



-hY.{p^ 



p; 



.(9) 



The Hamiltonian together with the constraints Eq. (||) 
can be viewed as some special Z^ gauge theory and is ana- 
lyzed below and in further detail in Appendices |^ and 
We find that generically this theory can have two decon- 
fined phases (in addition to the confined phase) with the 
phase diagram shown in Fig. ^j. Here we only describe 
the deconfined phase that obtains when all ring exchange 
couplings are large, i^ring S> ^. As explained below, this 
phase is a Z3 X Z3 deconfined phase. 

From here on, our focus is on the above Z^ ring ex- 
change Hamiltonian. We drop all tildes on Z^ fields ipr 
and superscripts on (which are now integers modulo 
3). Also, we consider a Zn generalizationo of the above 
Hamiltonian and carry out the analysis in the general 
case. This is done for clarity of notation. 

For Kj-ing ^ h, a, good caricature of the bulk ground 
state is given by the wavefunction 



\GS) = J2'\M), 



(10) 



where the primed sum is over all configurations {n^} that 
satisfy the constraints Eq. (^), i.e., J2reR^r = 0. 



3 



Ka 
h 





Fully deconfined 


Confined \ 


Partially*----..,^ 
deconfinedf Z^y7'~~~----........_^ 





C30 

K/h 



FIG. 2: Generic phase diagram of the Zm ring exchange on 
hexagonal lattice, Eq. (^). To explore the different deconfined 
states, we allow two different ring exchange couplings, Turing = 
Ka for A-type hexagons and ifring ~ K for B- and C-type 
hexagons (see Fig. 



Let us define flux through a given hexagon R 



(11) 



with the sign convention as in Fig. |^. The ground state 
has zero flux through each hexagon. Excitations above 
this ground state are Zjq vortices. For example, we can 
add a unit of flux through a given hexagon by applying 
a "string" operator as indicated in Fig. |[ The gap for a 
vortex carrying one unit of flux is 2_fi'ring[l — cos(27r/iV)]. 

Observe now (Fig. ||) that the hexagon centers R form 
a triangular lattice, which consists of three sublattices A, 
B, and C. Observe also that the flux-adding string oper- 
ator "steps" only through the same sublattice hexagons. 
We are thus led to the possibility of a topological dis- 
tinction between vortices on the different sublattices, in 
addition to the usual distinction between two vortices 
carrying different flux. Indeed, one can see that the topo- 
logically distinct situations can be characterized by say- 
ing that we have two species of Zjq vortices, say A and B 
vortices. Alternatively, if we want to preserve the sym- 
metry among the three sublattices, we can say that there 
are three types of vortices — A, B, and C vortices — 
but these are not independent and instead satisfy "fusion 
rules" such as 

($A = +l) X ($B=+l) ~ ($c = -l) ■ (12) 

This means that a nearby pair of -1-1 A and B vortices is 
indistinguishable from a — 1 C vortex (note also that the 
"states" on the left and on the right can be connected by 
local h terms in the Hamiltonian). 

Consider now introducing Z^ charges in the above 
gauge theory Eq. (|^), e.g., consider placing a pair of oppo- 
site ±1 charges on two hexagons i?i and R2: X^rei? ~ 
5rBi ~ This is appropriate when studying the 

charged sectors of the microscopic Hamiltonian Eq. (|l|) 
since the formal gauge structure represents the crucial 
coupling of chargons with the above Z3 degrees of free- 
dom. A charge can be added to a hexagon by applying 




FIG. 3: Hexagonal lattice on which the Zn ring exchange 
model Eq. (|) is defined. Nr = X^^g^^ir measures Zn "num- 
ber" on each hexagon; we can increase Nh on a given hexagon 
(★) by one by applying a "string" operator tp'^ iptl;^ tp . . . along 
the indicated path. $r = <j)i ~ <j)2 + 4>3 ~^ 4>4, + 4>5 — 4>6 mea- 
sures Zm "flux" through a given hexagon (123456); to fix the 
sign convention we always take 1 to be the lowermost hexagon 
site. We can increase the fiux through a given hexagon (•) 
by one unit by applying a string P'^ P'^ . . . along the 

indicated path. The lower-right corner of the figure shows the 
three sublattice structure of the lattice of honeycombs. 



a string operator as indicated in Fig. ^ From several 
perspectives, one can see that all such charges are de- 
confined in the -?Cring ^ h phase: Thus, in Appendix ^ 
we approach this "fully deconfined" phase starting from a 
"partially deconfined" phase, in which charges are decon- 
fined on one sublattice only. Also, this fully deconfined 
phase corresponds to the fully disordered phase in the 
dual global Z^ spin model of Appendix 

Similarly to vortices, we need to distinguish the charges 
on different sublattices. Again, as far as the gauge struc- 
ture is concerned, we have fusion rules such as 



(iV^ = +l) X (Nb = +i) - (iVc = -l) 



(13) 



Statistical interactions between the different particles 
are readily identified by studying the commutation prop- 
erties of the corresponding strings. These are summa- 
rized in Table | by specifying "gauge charges" of the dif- 
ferent Nji — +1 excitations with respect to the A, B, and 
C fluxes. Thus, iV^ = +l excitation carries gauge charges 
Qa — 0,Qb = — l,Qc = +lj i-e., it does not "see" A 
vortices, but when transported around a B or C vortex 
of unit strength, the wavefunction acquires an additional 



phase e 



-i2Tr/N 



or e 



i2ir/N 



correspondingly. 



This completes the particle description of the fully de- 
confined phase. The minimal description would be to say 
that we have A hexagon charges that see B hexagon Z^ 



4 



TABLE I: Gauge charges of the Nn = +1 excitations with 
respect to the A, B, and C Zn fluxes as defined by Eq. ( \L 1|) . 
Note that these are consistent with the fusion rules Eqs. ( |12| ) 
and 





Qa 


Qb 


Qc 







-1 


+1 




+1 





-1 




-1 


+ 1 






vortices and B hexagon charges that see A hexagon vor- 
tices. This is essentially the claimed Zjv x Zjv structure. 
Thus, we expect A^^-fold ground state degeneracy if the 
system is put on a cylinder, as can be verified by con- 
structing the corresponding ground states starting from 
the state Eq. (10) and threading vortices through the hole 
of the cylinder. 

It should be emphasized here that the above discussion 
assumed that the three-sublattice structure is respected 
by the boundary conditions. While it is clear that the 
bulk properties do not depend on this, there is an ad- 
ditional quirk when we consider topological degeneracy 
in a geometry that does not respect the three-sublattice 
structure. This is legitimate when all ring exchange cou- 
plings are equal. Consider, e.g., a cylindrical geometry 
with the circumference along the horizontal direction of 
Figs. |l| and ||. When an A-type particle is transported 
around the periodic direction, it does not return to its 
initial position, but rather becomes a B- or C-type parti- 
cle. It takes three turns for the particle to return to the 
original position. From table ||, such — +1 particle 
will not register any flux in this process. A detailed anal- 
ysis shows that for N mod 3^0 the ground state of the 
system in this geometry is nondegenerate. On the other 
hand, for TV mod 3 = the ground state is found to be 
three-fold degenerate, since in this case there is a com- 
posite object that returns to its initial state when trans- 
ported once around the cylinder and that senses some 
flux through the hole of the cylinder in the process. 

Returning to our microscopic bosonic model, the "par- 
ticle description" of the x phase is as follows: We 
have two species of Z3 vortices (with gap ~ if ring) and we 
have charged particles (with charge gap ~ U) that can 
be classified as carrying two distinct Z3 gauge charges, 
in addition to their fractional electrical charge. Finally, 
note that the Z3 x Z3 state is associated with the addi- 
tional symmetries in the hexagonal lattice ring exchange 
Hamiltonian but is protected by the same charge gap pro- 
jection, since any move within the uncharged state sector 
is necessarily a combination of hexagon ring exchanges. 



III. CONCLUSIONS 

We showed that it is possible to produce more com- 
plicated fractionalization patterns such as Z3 fraction- 



alization in relatively simple bosonic models. While 
the resulting fractionalized state turned out to be even 
more complicated than initially intended, the microscopic 
model was not too contrived. It is hoped that this work 
will encourage further searches for other exotic states. 
For example, can a non-Abelian fractionalized state be 
produced in a condensed matter system with a global 
symmetry only, short-range interactions, and in zero 
magnetic field? 
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APPENDIX A: PARTIALLY DECONFINED (Zn) 
PHASE 

To better appreciate the character of the deconfine- 
ment in the special gauge theory Eq. (^, we allow dif- 
ferent ring exchange couplings for different hexagons and 
consider particular parameter space with two such cou- 



plings: K^i 



Ka for the A hexagons and K-c[ 



Kc = K for the B and C hexagons. This is indicated 
schematically in Fig. || where the A hexagons are shaded. 
Note that by allowing the two couplings we implicitly as- 
sume that the boundary conditions on the lattice respect 
the three-sublattice structure; this is done throughout. 

We argue below that the ring exchange Hamiltonian 
has the phase diagram shown in Fig. with three phases: 
For h ^ K, Ka the system is in a fully confined phase. 
For K ^ h ^ ^/KaK the system is in a partially de- 
confined (Zn) phase. In this phase, the charges on the 
A hexagons are deconfined, while the charges on the B 
and C hexagons are confined. Finally, for K, Ka ^ h the 
system is in a fully deconfined {Zj^ x Zjv) phase with all 
charges deconfined. The phase diagram of Fig. ^ is also 
supported by the analysis of the dual global Zn spin 
model summarized in Appendix 

In what follows, we give a detailed description of the 
partially deconfined phase. As a representative of this 
phase, consider the Hamiltonian with Ka = 0, i.e., with 
ring exchanges around the B and C hexagons only (see 
Fig.^). In this case, there are additional conserved quan- 
tities: 



~ const (model with Ka = 0)(A1) 



for each hexagonal lattice link {rr') between two A 
hexagons A and A' (see Fig. This facilitates the 
analysis, since we can consider separately each subsec- 
tor specified by the corresponding eigenvalues {Caa'}- 
Note that the allowed {Caa' } are very much constrained 
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by the constraints Eq. (g) on the Ur themselves; however, 
we will not use the details of these explicitly. 

First of all, observe that the A hexagons in turn form 
a triangular lattice, while the links (rr') between such 
hexagons can also be viewed as the links of this "A- 
lattice", {AA') = {rr'). In a given subsector with fixed 
{CAA'}^ there remains one Z^r degree of freedom for each 
such link. It is convenient to work in the number basis 
and label these remaining link degrees of freedom by 



(0) ^_(„^,_„(o))^__^^,^^^ (A2) 



Ma-^A' —Ur — n. 



where is one particular instance: Caa' = '^r'^'' + 

(and our convention is that r € A and r' g A' — see 
Fig. H). Thus, AFaa' = ■M'a^A' are oriented fields on the 
links of the A-lattice. The subsector is now completely 
specified by the conditions 



J2 ^AA' 

A'eA 







(A3) 



which are the neutrality constraints Eq. (|g) for the A 
hexagons. 

The action of the Hamiltonian Eq. (^ in this subsec- 
tor is readily described in terms of the new variables. 
Thus, the transverse field (h) terms are diagonal in the 
new number variables, while the B and C hexagon ring 
exchanges simultaneously raise (or lower) the three ori- 
ented number fields circulating around the corresponding 
A-lattice triangular plackets. Writing the raising opera- 
tor for a given link number variable Maa' as e''^^^' , the 
resulting Hamiltonian is 

H[C] ^ -KY^ ^ j^^^^^ (-^4^ 

A 

J2 (r^A'e-*(2-/^)-^--'+h.c.) , 



(AA') 



where 



Together with the constraints Eq. (A3), this is precisely 
the conventional Z]\j lattice gauge theory defined on the 
triangular A-lattice but with link-dependent Taa' specific 
for the particular subsector {Caa'}- We can now use the 
conventional wisdom to characterize each such subsector 
and in turn the full hexagonal ring exchange Hamiltonian 
with Ka = 0. 

When /i = 0, all the different subsectors are degener- 
ate. The lowest energy state in each such subsector has 
the energy of —2K per triangle and is an equal weight 
superposition of all possible con figurations of Maa' that 
satisfy the constraints Eq. ([A^). Nonzero h eliminates 
this degeneracy and selects one particular subsector, 
namely with all Caa' = 0, as containing the true ground 
state of the full Hamiltonian with Ka = 0. Indeed, treat- 



ing V AA' perturbatively, the lowest energy in a given sub- 
sector is 



ITaa' 



A 



^^^^^ 2K[l - cos(2^/7V)] 



(A6) 



where for simplicity we assumed that the system has no 
boundaries. It is now clear that for small nonzero h the 
ground state of the full ring exchange Hamiltonian with 
Ka = is in the subsector with all Caa' — 0. The sub- 
sectors that are closest in energy have the smallest num- 
ber of nonzero Caa' and can be characterized as having 
alternating CaqA' = +1 and CaqA' = — 1 values on the 
six links to a given hexagon Aq (this subsector is obtained 
from the ground state subsector by applying the hexagon 
ring exchange around the hexagon Aq). The energy gap 
to these subsectors is /K. We see that we have a 
peculiar situation where a nonzero transverse field h is 
needed to stabilize this Z^q deconfined ground state; this 
is because we are competing here against the Zjq x Ziq 
deconfined state that is obtained for large Ka,K. 

We are all set to discuss confinement of charges in the 
model with Ka = 0. The above analysis was carried 
out in the uncharged sector but is readily extended to 
the charged sectors. First, consider placing a pair of op- 
posite charges on two B hexagons: X^rGS ~ ^'^^ 
= — 1. Proceeding exactly as before, one is led 



to consider different subsectors (of this charged sector) 
specified by {Caa'}- In each such sub sector, the Hamil- 
tonian has precisely the form Eq. (A4) with the number 
variables satisfying precisely the constraints Eq. (|A3| ). 
All information about the two charges is encoded in the 
allowed configurations {£aa'}, and one can clearly see 
that Caa' ^ at least on a string of A-lattice bonds 
connecting Bi and -B2. From the earlier arguments, the 
energy cost of introducing two such charges is then pro- 
portional to the length of this string, i.e., such charges 
are confined with the string tension h? / K . 

Consider now placing a pair of opposite charges on two 



A hexagons Ai and A2: X] 



rGAi 



-1 and X 



r&A2 



— 1. The analysis of the subsectors {Caa'} will be some- 
what different in this case. For each such subsector in this 
charged sector there corresponds a subsector in the un- 
charged sector having exactly the same {Caa'}- It is con- 
venient to "measure" each charged subsector relative to 
the corresponding uncharged subsector. This is achieved 
by defining link variables Maa' via Eq. ( [A2| ) using an un- 
charged instance {rir*^^} of {Caa'} (i-c, 



,(0) 



for each R and Laa' = ^^T■°^ + n^V'' for each {AA')). 
In ea ch subsector, the Hamiltonian again has the form 
Eq. (A4) when written in these link variables, which 



,(0) 



reR. ' 







now satisfy new constraints ^ 



A'eA- 



A'eA^■^A^A 



+1 and 



A2A' 



= — f. This corresponding precisely to 
introducing two charges in the corresponding A-lattice 
gauge theory. Clearly, for large enough K ^ h, these 
charges will be deconfined. 

We now have essentially complete description of the 
partially deconfined phase. Thus, one can readily iden- 
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tify the Zn vortex excitations of the A-lattice gauge the- 
ory with Zj^ vortices on the B and C hexagons. These 
vortices will have usual statistical interactions with the 
deconfined charges on the A hexagons. Also, as should 
become clear by reviewing the above discussion, we can 
essentially account for the different subsectors {£aA' } by 
saying that there are additional particle excitations living 
on the A hexagons obtained from the ground state by the 
action of the corresponding A hexagon ring exchanges. 
These new particles have a "mass" of &h?/K and have 
no statistical interaction with the other particles. 

We can now consider what happens when we allow 
nonzero Ka- As discussed above, the ring exchanges 
around the A hexagons introduce mixing between the 
different subsectors. However, as long as Ka is much 
smaller than the corresponding gap ~ h? /K, the par- 
tially deconfined phase survives and is characterized by 
the same particle description. 

Once Ka is sufficiently large, the system enters the 
fully deconfined phase described in the main text. 

APPENDIX B: DUAL GLOBAL Zn SPIN MODEL 

Here we summarize dual perspective on the hexago- 
nal lattice ring exchange Hamiltonian Eq. (^). We work 
directly in the Hamiltonian language. Simple counting 
shows that the dimensionality of the physical Hilbert 
space is consistent with having one Zj^ degree of free- 
dom per hexagon. Let us define 

= ^\^2i^lMli'6 , T+ EE (T^)t , (Bl) 

where we use the same sign convention as in Fig. |^. Let 
us also define 

Vl^^\{P+ , (B2) 

where the product is along the vertical path that reaches 
R as in Fig. ||. Note that the path "steps" through 



the same sublattice hexagons. [If we were to take some 
other such path, we would need to replace some P+ with 
. The total product is path-independent due to con- 
straints Eq. (||).] Vjj can be thought of as a vortex cre- 
ation operator. 

We now interpret as a Zj^ spin variable. It is easy 
to verify that is the corresponding conjugate variable 
(i.e., raising operatorEJ): 

VIt+ = e^'^'^'^T+vl . (B3) 
Also, we can readily "solve" for P+: 

The dual Hamiltonian is 

H = -KJ2{T^ + h.c.) ^hY^iV^Vlvl + h.c.)(J35) 

R A 

which is a global Zjv spin model with three-spin interac- 
tions. A little thought shows that the model has in fact a 
Zn X Zn global symmetry corresponding to independent 
global rotations of the spins on two of the three sub- 
lattices. Note also that the three-spin interaction around 
triangles promotes ordering of the spins on the same sub- 
lattice. This is because two neighboring sites A and A' 
on the same sublattice share a BC side in the respective 
triangle interactions AABC and AA'BC. 

The global model clearly has a fully disordered phase 
for K ^ h. In the original ring exchange Hamiltonian, 
this corresponds to all vortices being gapped, and we 
obtain the Zn x Zn fully deconfined phase. Varying the 
A hexagon ring exchange coupling Ka independently, for 
sufficiently small Ka and large K the system can clearly 
order on the A sublattice (i.e., A vortices condense), but 
remain disordered on the B and C sublattices. This is our 
partially deconfined Zn phase. 
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